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ABSTRACT 



The structure of a hydro dynamic, double transonic flow driven by neutrino 
annihilation in the polar region of a Schwarzschild black hole is computed for 
different energy deposition profiles. The requirement that both, the inflow into 
the black hole and the outflow to infinity pass smoothly through their sonic 
points fixes the stagnation radius and stagnation pressure. The asymptotic power 
of the outflow is shown to be the integral of the energy deposition rate above 
the stagnation radius. The outflow production efficiency depends on the energy 
deposition profile, and is generally higher for shallower profiles. Using recent 
calculations of the neutrino annihilation rate, we estimate that over 50 percents 
of the total energy deposited above the horizon can emerge in the form of a 
relativistic outflow at infinity. The continuous creation of plasma during the 
expansion of the outflow leads to generation of a large specific entropy. This has 
important implications for the prompt photospheric emission in GRBs. 



1. Introduction 

The relativistic outflows producing the prompt and afterglow emissions in GRBs are 
commonly thought to be powered by hyper accreting black holes. The high Lorentz fac- 
tors inferred from energy considerations, compactness arguments and afterglow models, 
r ~ 10^ — 10^, require extremely low baryon load at the outflow injection point, which 
seems difficult to achieve in disk outflows (Levinson 2006; Barzilay & Levinson 2008; Met- 
zger et al. 2008). A natural way to avoid baryonic contamination is to accelerate the outflow 
in the region above the horizon of the black hole (Levinson & Eichler, 1993). Two popular jet 
production mechanisms, that are widely discussed in the literature, are magnetic extraction 
of the rotational energy of a Kerr black hole (Blandford & Znajek 1977), and pair creation 
on horizon threading field lines by annihilation of neutrinos that emanate from the accretion 
disk surrounding the black hole (Eichler et al. 1989, Popham et al. 1999; Asano & Fukuyama 



^School of Physics & Astronomy, Tel Aviv University, Tel Aviv 69978, Israel 



-2- 



2001; Birkl et al. 2007; Zalamea and Beloborodov 2011; hereafter ZBll). The former mecha- 
nism requires nearly maximal rotation of the black hole and sufficiently high magnetization. 
The latter mechanism can operate also in a Schwarzschild spacetime, although extremely 
high accretion rates ( ^ 3Mq s~^) are required in this case to account for the characteristic 
GRB luminosities (ZBll). It is worth noting that the presence of a rapidly spinning Kerr 
black hole may not be sufficient to guarantee effective production of a Blandford-Znajek 
(hereafter BZ) outflow; overloading of magnetic field lines by relativistic pairs may lead to a 
complete shutdown of this mechanism, in which case one has to rely on the neutrino-driven 
outflow to power a GRB. However, the accretion rate required to power a typical GRB by 
the latter mechanism can be considerably reduced when the specific angular momentum of 
the black hole, a, is high enough. The reason is that the neutrino luminosity emitted by the 
disk and, hence, the rate at which energy is deposited above the horizon by vv annihilation, 
increases sharply with increasing a (ZBll). 

A question of interest is what fraction of the total energy deposited above the horizon 
via vv annihilation emerges at infinity. To address this question, we constructed a model for 
the double-transonic fiow established in the polar region, that incorporates energy injection 
by the external neutrino source in a self-consistent manner. 

2. Model description 

Neutrino annihilation in the baryon free region above the horizon generates a relativis- 
tically hot fiuid consisting of a mixture of e^ pairs and radiation in equilibrium. Part of 
this fiuid is accelerated outwards by pressure gradient forces, and the other part falls into 
the black hole, owing to the strong gravitational force exerted on it. In what follows, we 
consider the structure of the double transonic fiow thereby produced. Schematic illustration 
of the model is shown in figure [H We restrict our analysis to the Schwarzschild spacetime, 
and suppose that the fiow is stationary, radial, unmagnetized and non-rotating. For clar- 
ity we consider baryon free plasma, although our analysis applies also when baryons are 
present, provided the enthalpy per baryon is sufficiently large {w/rib » 1). The idea is to 
find a self-consistent solution that starts from a stagnation point, where the fiow velocity 
vanishes, and crosses two sonic points; an inner one located in the infiow section below the 
stagnation point, and an outer one located in the outfiow section above the stagnation point. 
As shown below, for a given choice of energy deposition profile the stagnation radius and 
stagnation pressure are uniquely determined by the requirement that the solution passes 
smoothly through both sonic points. 
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2.1. Energy deposition rate 

The annihilation of neutrinos emitted from a hyper-accretion disk was calculated in a 
number of recent works, under different simplifying assumptions (e.g., Popham et al. 1999; 
Asano & Fukuyama 2001; Birkl et al. 2007; ZBll). The most detailed analysis is presented 
in ZBll, who adopted the relativistic disk model of Chen and Beloborodov (2007) as the 
neutrino source, and employed a geodesic-tracing method in Kerr spacetime to evaluate the 
local energy-momentum deposition rates by the reaction vD — > e''"e~ in the vicinity of the 
black hole. 

The local energy-momentum deposition rates, Q%, were computed in ZBll in the frame 
of a zero angular momentum observer (ZAMO). Quite generally, these rates are functions of 
r and 9, when expressed in Boyer-Lindquist coordinates, = {t,r,9,(f)). Since we restrict 
our analysis to a conical flow, we shall adopt a local energy deposition rate of the form: 

QUr) = Qofir/rn), (1) 

where Qi^ir) represents the angle- averaged rate at radius r, th is the horizon scale, and 
/(I) = 1. From figures 2 and 3 in ZBll we estimate f{x) ~ for a = 0.95 and 

f[x) ^ x^^'^ for a = 0, in the range < r < 30m delineated in the figures, where m is 
the black hole mass in geometrical units and a its specific angular momentum. The profile 
should steepen as the radius increases, approaching f{x) oc x^^ at radii much larger than the 
size Ri, of the neutrino source (e.g., Goodman et al. 1987, Qian & Woosley 1996). For the 
disk model outhned in Chen & Beloborodov (2007), Ri, can be identified with the ignition 
radius of the disk, rjp„, within which neutrino emission is efficient. 

The net energy deposition rate, as measured by a distant observer, is given in terms of 
the ZAMO rates Q% by (ZBll) 

= / f ^k= - -^Qt) V^drd9dcl>, (2) 

Jr>rH \ V9H J 

here g^i, are the metric components of the Kerr spacetime in Boyer-Lindquist coordinates. 
The results exhibited in figure (4) of ZBll (see also their Equation (22) for a fitting formula) 
indicate that E^^p depends sensitively on the accretion rate m (given in units of Mq/s) and 
the spin of the black hole a, with E^p ~ 5 x lO^^m^^^ erg s~^at a = and E^p ~ lO^^m^^"^ 
erg s~^ at a = 0.95, for a black hole mass Mbh = ^Mq, and accretion rates in the range 
Mign < M < Mfrap- At M < Mign ucutrino emission is severely suppressed by virtue of the 
low disk temperature. At Mtrap < M the neutrinos are trapped in the disk and advected 
into the black hole, whereby the neutrino luminosity saturates. 
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2.2. Flow equations 

The stress-energy tensor of a purely hydro dynamic flow takes the form 

T"^ =^/;u"u^+p^"^, (3) 

here m° is the four-velocity measured in units of c, w is the specific enthalpy, p the pres- 
sure, and g^jj = diag(— a^, r^, sin^ 6'), = 1 — 2m/r, is the metric tensor of the 
Schwarzschild spacetime. 

The dynamics of the flow is governed by the energy-momentum equations: 

^ ' ^T"''),« + rVT^^ = /, (4) 



V-9 - - 

where denotes the source terms associated with energy-momentum transfer by the external 
agent, and F^^, denotes the affine connection. By contracting gi3^ with Equation (jl]), using the 
relation {\/—gg'^^),a + y/—g^^,yg^'^ = 0, taking the t component and noting that T ^tvU^u'^ = 
for a stationary flow, one obtains 

^wu°'ut),a = qt- (5) 



Likewise, contracting with Equation (jl]), using the identity V^^^^upu^u'' = —Ui3u''{u'^)^^, 
yields 

^WU°')^a - u"'P,a = -U^q" ■ (6) 



Combining Equation (jG]) with the thermodynamic identity d[wV^ — Vdp = kTd{Vs), where 
V is the volume of a fluid element, T its temperature and and s its specific entropy, yields 
the change in s: 

kT 

^su°')^a = -u^q". (7) 



Henceforth, we consider a conical flow, for which = (m*,u^,0,0), m^(9^ = u'^dr, and 
denote by Up = u^/a the poloidal velocity, and by 7 = (1 -|- UpY^"^ and v = Up/'j the 
corresponding Lorentz factor and 3-velocity. The convention here is that < on inflow 
lines and -Up > on outflow lines. Since the created pairs are relativistic, we adopt the 
equation of state w = Ap. Furthermore, we neglect momentum transfer to the flow by the 
created pairs, so that q'^ = (g*, 0,0,0). The source term g* is measured in the frame of a 
distant observer, and can be expressed in terms of the local energy deposition rate adopted 
in Equation ([1]) as eg* = V~9^*Qiu — Qof{^)/(^- Under the above simplifications. Equations 
(ISl) and (16]) reduce to 



jUpAp (1 — X ^y/'^ X x{x — 1) 
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and 

dxlnp= ^ — - {l + v'^)d:,\iaup, (9) 

7Mp4j9(l — X ^)^'^ X x{x — 1) 

in terms of the dimensionless radius x = r/r^, = 2m, and the normahzed pressure p = 
p/{Qotd), td = Vg/c being the dynamical time of the black hole. Equation ([8]) has critical 
points at f = ±Cs, where Cg = is the sound speed. 

At the stagnation point x = Xgt, where Up = 0, the above equations yield 

= - ^ (11) 

where pst = p{xst) is the stagnation pressure. Thus, for a given choice of /(x) the solution 
is fully determined once Xgt and Pst are known. 

The regularity conditions at the sonic points, obtained from Equation ([H]), read 

VSxlil - x',^f'^f{x,{) = -2(2x,i - 3)p,i, (12) 

and 

V^xUl - x-^Y'^f{x,2) = 2(2x,2 - 3)p,2, (13) 
denoting the sonic point of the inflow (outflow) by Xci(xc2), and noting that Up = —l/\/2 at 
Xci, and Up = l/\/2 at Xc2- Evidently, Xd < 3/2 and Xc2 > 3/2. The existence of two sonic 
points is a consequence of energy injection. When f{x) = the system has only one critical 
point at Xc = 3/2. 



2.3. Asymptotic power and entropy of the outflow 

Well above the sonic radius, at x >> Xc2, the power of the outflow is given by Lj^o = 
Itic JII° T^^r'^ sin 6 dO = 27r(l — cos6'o)w7Mp, here 6q is the opening angle of the flow. Inte- 
grating Equation from the stagnation point Xgt to infinity, using the boundary condition 
Up{xst) = and noting that qt = gtt(f = — a^g*, one obtains 

/»00 /"OO 

Ljoo = 27r(l-cos^o) / cqtr'^dr = 2n{l- COS do)Qorl i {I - x'^f/^ f{x)x^dx. (14) 

Let us define the outflow production efficiency e to be the fraction of Eyp that emerges at 
infinity. From Equation ([2]) with gtif, = 0, y/—g^ = a~^, and Equation f|T^ . one finds 

_L,^_ jZ{i-x-'y^'f{^Wdx 

' JZ{i-x-^y/^f{x)x^dx' ^ ' 



- 6 - 



The rate at which entropy is carried by the flow at some radius r is obtained by inte- 
grating Equation ([7]) from the stagnation radius r^^ to r > r^j: 



= 27r(l - co^e^ycu's = 27r(l - cos^o)^^^ T ip~^'^ f{x)x'^dx\ (16) 

where Tq is a fiducial temperature defined bjfl To = {UQotd/llay^*. If the outflow also 
carries baryons at a rate A*";, = 27r(l — cos 6o)nhU^r'^ , for a baryon density n?,, then the 
dimensionless entropy per baryon can be expressed as 

= if = ^^^F{x), (17) 
Nb at kTo 

with = Ljoo/ {NbTripC^) being the maximum Lorentz factor of the outflow, and 

Ix 'lP~^^^f{x')x''^dx' 



F{x) 



(l-x'-i)V2/(a;')x'W 



3. Results 

Equations (IE]) and ([9]) have been integrated numerically using energy deposition profile 
of the form /(x) = x~^ . The pressure at the inner sonic point, Pci, is used as a free parameter 
that we adjust to converge to the desired solution. The integration starts at the inner sonic 
point Xci , where Up{xci) = — 1 / \/2, and repeated iteratively by changing the value of Pci , until 
a smooth transition across the outer sonic point is obtained. The value of Xd is computed, 
in each run, from the regularity condition (fT2l) . The stagnation radius Xgt is then found 
numerically from the condition Up{xst) = 0. 

Solutions are exhibited in figure [2] for different values of b. The horizontal dashed lines in 
the left panel mark the sonic velocities of the inflow {up = —1/ a/2) and outflow {up = 1/ a/2). 
As seen, in all cases the sonic point of the outflow is located well above that of an adiabatic 
flow (xc = 3/2). This is a consequence of the injection of energy over extended scales, that 
loads the flow and delays its acceleration. The values of the stagnation radii found from 
the integration are Xgt = 1.57, 1.70, 1.81 and 1.95 for b = 5,4,3.5, and 3, respectively. For a 
comparison, note that the radius at which the escape velocity, Vesc = \^rjr = equals 
the sound speed, Cs = l/VS, is x = 3, slightly above the stagnation point. The corresponding 



^The pressure of a relativistically hot plasma at a temperature T, consisting of electrons, positrons and 
radiation in equilibrium, satisfies p = llar*/12, from which we obtain T = Tqp^^^. 
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efficiencies, computed using Equation f lTSj) . are e = 0.58, 0.73, 0.83 for b = 5,4, 3.5. For 6 = 3 
the energy deposition rate diverges logarithmically and needs to be cut off at some radius 
x,^. We arbitrarily invoked Xi, = 100, whereby e = 0.93. 

By employing Equations ( IT4|) and (fT7|) . we can express the entropy per baryon in the 
form: 

^ 1.3 X 10=,, (^)(^)-"(^||) (19) 

where Lj^o = 2Ljoo/(l — cos 6*0) is the isotropic equivalent luminosity, Mbh is the mass of the 
black hole, and rji, = 1.06,0.84,0.62 for b = 3.5,4,5, respectively. A plot of F{x) is shown 
in figure [3l The asymptotic value of F{x) represents enhancement of the entropy relative to 
that produced in an adiabatic outflow injected from a radius Rq ~ r^. 

3.1. Asymptotic behavior 

As pointed out in Section 12.11 sufficiently far out the energy injection profile must 
steepen, approaching f{x) oc x~^ at x » Ru/^m, where Ry is roughly the size of the 
neutrino source (i.e., the disk radius within which neutrinos are emitted). The first term on 
the right hand side of Equation ([8]) then becomes negligibly small, and the outfiow enters 
the adiabatic regime, whereby Equation (|8]) can be solved analytically in the limit f ~ 1 and 
/(x) = 0. The solution reads: Up{x) oc x'^/'^ /{x — lY^"^, and reduces to the well known result 
Up (X X, p (X x~^, at a; > Xc2 » 1- From Equations dH]) and it is clear that the solution 
will approach this asymptotic behavior once 6 > 4, as indeed confirmed in figure O 

Near the horizon we obtain the analytic solution Up oc — l/(x — 1)^/^ for any value of 
b, in accord with the numerical solutions exhibited in figure [2] , indicating that the infiow 
moves along radial geodesies as it approaches the horizon. This is a consequence of the fact 
that near the horizon the dynamics of the infiow is dictated by the gravitational attraction 
of the black hole. 



4. Conclusions 

We have computed the structure of a double transonic fiow generated above the horizon 
of a Schwarzschild black hole by annihilation of neutrinos emitted from a hyperaccretion disk. 
We have shown that for a given choice of the energy deposition profile, there exists a unique 
solution that passes smoothly through the inner and outer sonic points. The stagnation point 
was found to be located slightly below the escape radius of sonic material - in the range 1.5 - 



8 



2 Schwarzschild radii, depending on the energy deposition profile. The asymptotic power of 
the outflow is given as the integral of the energy deposition rate above the stagnation radius. 
The outflow production efficiency was found to be typically large, with over one half of the 
injected energy emerging at infinity. 

The continuous injection of energy pushes out the outer sonic point and delays the ac- 
celeration of the outflow above the sonic point. The linear acceleration phase commences 
once the energy deposition rate becomes small and the outflow enters the adiabatic regime. 
The radius at which this happens depends on the energy deposition profile, and may exceed 
100m at large accretion rates. The deposition of energy over extended scales leads to gener- 
ation of a relatively large specific entropy. As a consequence, in models wherein the prompt 
emission originates from the photosphere (e.g., Ryde & Peer 2009; Lazzati et al. 2009; Peer 
et al. 2012; Levinson 2012; Beloborodov 2013), the thermal peak may be located at an 
energy lower than previously estimated. For instance, if the prompt emission is produced by 
sub-photospheric, radiation mediated shocks (Bromberg et al. 2011; Levinson 2012), than 
the observed temperature behind the shock is 



where Ug is the shock 4- velocity (i.e., the 4- velocity of the upstream fluid, as measured in 
the shock frame), and F^o is the asymptotic value of the function F{x) exhibited in figure O 
Consequently, the peak of the spectrum produced by mildly relativistic shocks (m^ ~ 1) can 
be located at an energy Epeak ~ ^kTobs ~ 200 keV, significantly lower than those exhibited 
in Levinson (2012). Additional dissipation and entropy generation may result from the 
interaction of the outflow with the stellar envelope, via formation of re-collimation shocks 
(Bromberg & Levinson 2007). 

In situations where the central black hole is rapidly rotating, the outflow can be produced 
by the Blandford-Znajek mechanism. The resultant Poynting power is expected to exceed 
the net power deposited above the horizon by the reaction vv — )■ e^e~ (e.g., Kawanaka et 
al. 2012). However, activation of the BZ mechanism requires sufficiently high magnetization 
near the horizon, and it is therefore anticipated that overloading of magnetic field lines by 
relativistic pairs may lead to a complete shutdown of this process. In the latter case, the 
outflow will be driven by the pressure of the injected plasma, as described in Section^ rather 
than by magnetic extraction of the spin energy of the black hole. Interestingly, relatively low 
accretion rates may be favorable for producing powerful outflows. The critical load above 
which activation of the BZ process is prevented is a key issue, currently under investigation 
(Globus & Levinson, in preparation). 




keV, 
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Fig. 1. — Schematic illustration of the double-transonic flow model. 
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Fig. 2. — Velocity (left panel) and pressure (right panel) profiles for f{x) = {x = r/rg) 
and different values of b, as indicated. The region above (below) the dotted line = in 
the left panel corresponds to the outflow (inflow). The horizontal dashed lines delineate the 
sonic velocities of the inflow (lower line) and outflow (upper line). 
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Fig. 3. — A plot of the function F{x) defined in Equation ( IT8l) . for f{x) = x ^ and different 
values of b, as indicated. 



